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ABSTRACT 

The determination of the optimal distribution of aimpoints is 
examined for weapons that fire fragmenting projectiles against mobile 
targets. A finite difference approximation which reduces the problem 
to a mathematical programming problem is developed. Computational 
considerations for this nonlinear programming problem are discussed. 
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I. NATURE OF THE PROBLEM 



Consider a situation where an artillery battery or naval gunfire 
support ship desires to fire at a moving land target. If the original 
position, speed, and direction of the target were known, the guns 
cotild be aimed so that the projectiles would impact at a calculated 
future location of the target. This future location is fairly easily 
calculated in this simple situation. The calculations become more 
difficult whenever any of the factors of original position, speed and 
direction are unknown. In these cases a probabilistic determination 
of future target location becomes necessary. It is then possible to 
state only that the target is more or less likely, in a probabilistic 
sense, to be at a certain location than at another. 

This paper will examine the problem of determining the optimal 
distribution of aimpoints for a weapon which fires a projectile against 
a moving target whose speed is known but whose original position and 
the direction of movement are not. 

A model for the location of a target under the above circmstances 
was developed by B. Koopman in 194^ (Ref. l), and is described herein. 
AlthoTigh formulated in a Naval setting, Koopman 's model equally well 
applies to a moving target that has been detected by a for\^ard observer 
the exact position of the target being unknown. All that is known is 
that the target is more likely to be at a point 0 than at any other 
point. The target may not be "at 0, however, but only within a short 
distance of 0, all points the same distance r from 0 being equally 
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likely. The probability that the target is in an area dA at a distance 
r from the target is defined as P(r)dA. Koopman assumed that the 
situation could be approximated by the circular normal distribution, 

P(r) = exp -[-rVss^j , (1) 

2rrs 

2 

where s is the variance in the original target location. 

The speed of the target, u, is assumed to be known but the direction 
of the target movement is unknown, all directions being equally likely. 
After t units of time the situation may be as shown in Figure 1 . 



iA (NiriACUY 




Figure 1 . Entry of target into dA 

Koopman determined the distribution of moving targets about a point 0 
to be 

P(r,t) = exp |-(r^+u^t^)/2s^j (^) (2) 

2tts s 

where I is the modified Bessel function of the first kind and zero 
. o . 

order. The graph of P(r,t) for different t is shown in Figure 2. 
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Observe that the probability spreads outward with time so that the 
target is most likely to be in an expanding ring about 0. 




Figure 2 . The distribution of moving targets about a point 0. 
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II. MTHEMATICAL MODELS OF WEAPONS FIRING 



FRAGMENTING PROJECTILES 

It was necessary to examine the mathematical formulation of a 
weapon that fires fragmenting projectiles in order to determine the 
distribution of aimpoints at a moving target. The basic models for 
such weapons usually include four probability distributions. 

The probability that a target located at coordinates x^, y^ is 
killed by a round impacting at coordinates x,y is defined as P^. (x, -x,y -y) . 
This probability may be considered as the lethality function. The 
lethality function is the conditional probability that a target at x , 

y 

y^, is killed given that a round impacts at x,y. It may take several 
forms of which exponential, linear and "cookie-cutter” are familiar. 
Reference 2 discusses lethality functions. 

The distribution of impact points about the point of aim is defined 
as P^(x,y). This distribution is caused by meterological effects on 
the projectile and ballistic effects inherent in the weapons system. 

Some models have incorporated the distribution of impact points about 
the point of aim into the lethality function. It is then possible to 
define the probability of killing a target located at coordinates x^, y^ 
with a round aimed at coordinates x , y , This probability is defined 

3 / Si 

The probability distribution of aimpoints for a target is dependent 

upon the nature of the target, its movement and location, and the type 

of weapons system employed. It is defined as P.(x ,y ), 

A 3, 3/ 
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The probability distribution of the location of the target is 
defined as P^(x^,y^). 

These distributions have been combined to model weapons systems. 

This paper examines three of these models to illustrate the modeling 
techniques and to provide backgroimd for examining the problem of 
determining the distribution of aimpoints of a weapon firing fragmenting 
projectiles against a moving target. Although a different notation is 
used in each of the original source documents, a standard notation has 
been adopted for purposes of presentation in this thesis. 

A. THE GROVES MODEL 

This model was developed by Groves in Ref. It is a simple model 
that does not consider a distribution of aimpoint; i.e. the aimpoint 
is fixed. If P,. (x. -x,y -y) is the lethality function representing the 
probability that a target located at coordinates x^, y^ is killed by a 
round impacting at coordinates x,y, and P^(x,y) is the distribution of 
the impact points; then the probability of killing a target at x^ , y^ 
by a round aimed at coordinates x , y is 

Si Si 

00 00 

^K^V^a’^f^a^ =11 Pj^(x^-x,y^-y)p^(x,y) dx dy. (3) 

_00 vOO 

The probability that a target survives one round is then l-P„(x,-x ,y,-y ), 

0 Si o Si 

and the probability of surviving N rounds all aimed at the same aimpoint 

IT 

is (l-P^(x.-x ,y.-y )) . The probability that the target is killed by 

K. tJ 3/ t/ 3, 

, , . N 

any of the N rounds is l-(l-P^^(x -x ,y,-y )) . This expression can then 

iv *tJ 3. *U 3/ 

be integrated over the area of the target, T, to obtain the expected 
fractional kill, K 

it = Cl-(l-Pjj(x^-x^,yj-y^))’')&^ dyj. (4) 

T 
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The substitution of equation (5) into (4) yields 

CO 00 

•p «.oo «.oo 

Groves further reduces this expression to one that is more suitable for 
hand-computing the expected fractional kill. It is noted that the 
simplicity of this model is a result of the fixed aimpoint as seen by 
comparing the Groves model with the Weiss and the Breaux-Mohler models. 

B. THE BREATJX-MOHLER MODEL 

This model was published by Breaux and Mohler in Eef. 4 * The model 
contains all of those distributions discussed earlier. It is used by 
Breaux and Mohler to compute the expected fraction of a target damaged 
by a salvo of fragmenting projectiles all aimed at the same aimpoint. 

They define the probability that a target located at coordinates 

X,, y, is damaged by a round impacting at coordinates x,y as P^. (x, -x,y , -y) , 

t/ X L X X 

Pj(x,y) is the density function describing the distribution of impact 
points x,y about the aimpoint x , y . It is assumed that all N rounds 

Si Si 

of a salvo are aimed at the same aimpoint. The probability a target 
survives N rounds is 

(1- J J Pj^(x^-x,y^-y)Pj(x,y) dx dy)^, . (6) 

^00 —00 

and the probability of damage over all impact points is 
00 00 ’ 

1 - (1- J J Pj^(x^-x,y^-y)Pj(x,y)dx dy)^. (?) 

—00 00 

Breaux and Mohler assume that the target is distributed over the target 

area, T, as Pm(x, ,y.), and the aimpoint itself is a random variable 

Jl X X 
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distributed as The expected fraction of the target damaged 

is determined to he 

f = JJ* J ! i' - i* <3y) 

^ ^00 ^00 m—00 _co 

This expression is reduced using a binomial expansion to, 

f . S /J / J [ J J P (x -x,y -y)P (x,y)ax^]^ 

Breavix and Mohler further reduce eqoiation (8) using Jacohi polynomials 
to produce an expression that can be used for determining the expected, 
fractional kill. 

C. THE WEISS MODEL 

This model was developed by Weiss in Ref, 5* It will be discussed 
more thoroughly than the Groves and Breaux-Mohler models because the 
model is used later in the paper, Weiss assumes a target of n men that 
is distributed as P^(x^,y^). The probability that a man is in the 
small area dx^ dy^ is P^(x^,y^)dx^dy^ and 

/ / ho) 

..00 ^00 

The expected number of targets in a small area is nP,j|(x^,y^)dx,|.dy_j., 

The probability that a round aimed at the aimpoint x , y will kill a 

Si Si 

target is P,,(x,-x ,y,-y ). Weiss assumes that N roiinds are fired at 
K t a t a 

the target and the ith round has an aimpoint x . » y • . The probability 
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that a target at x^, siuTvives all N rounds is 



N • - • 

1=1 



and the probability that a target is killed is 



( 11 ) 



( 12 ) 



The probability that there is a target at x^, y^ and that is killed is 
The expected number killed (K) is 

CO CO 

K = n J J P^(x^,y^)(l-q(x^-x^,y^-y^)jdx^dy.^, 

^00 _co 

00 00 00 00 

«00 00 ”00 _00 

00 00 

= ^ [l- J S PT^^t’^t^^^V^a’^t-ya^'^t'^^tJ* 



^00 ^00 



H“K 

The expected fraction surviving $ is , or 



00 00 



* - I - S [l- I I PT(^f5't’l(^t-’'a'^t-ya) 

^00 .00 

00 00 

- J" J '’T(*t’5't>l(v^a>yt-5'a5 



( 15 ) 



.00 .00 



If each round has the same aimpoint, then 



N 






( 16 ) 



.00 .00 



and iising a binomial expansion 

.00 .00 i =0 



( 17 ) 



.00 .00 



11 



Asaxuning that the aimpoints themselves are distributed as P,(x ,y ), 

3i 3* 

the expected fraction of targets surviving is 



00 00 
»00 mOO 



( 18 ) 



N 



00 00 






^00 ^00 



0=0 



^00 «»00 



' dx^dy^ dx^dy^. (I9) 

D. COMPARISON OP BREATJX-MOHLER AND WEISS MODELS 

It is easily shovm that the Breaux-Mohler and Weiss models give the 
same results. Substituting equation (5) into equation (s) gives the 

expected fraction damaged as 
00 00 

IP —00 — C30 

Substituting equation (l6) into equation (is) gives the expected fraction 
of targets surviving as 

00 00 00 00 

^=11 V^a’^a) I I ‘^^t'^yt'^a'^^a^^^) 

_00 _00 mOO mOO 

It is clear that these models are simply complements of each other. 



III. THE DETERMINATION OP NECESSARY AM) 

SUPPICIENT CONDITIONS 

The mathematical models described in chapter II can be used to 
determine necessary and sufficient conditions for determining the 
aimpoint distribution that maximizes the expected fraction of the 
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target destroyed or minimizes the expected fraction of survivors, 
depending upon the model used. Reference 6 discusses mathematical models 
of hit probabilities and the techniques of maximizing salvo kill 
probabilities. Weiss (Ref. 7) attacks this problem using techniques of 
Svesnikov (Ref. 8) and Morse and Kimball (Ref. 9)* His approach to the 
problem of maximizing the salvo kill probability (or minimizing the 
fraction of survivors) required the application of the calculus of 
variations. This application of the calculus of variations can be 
traders tood by examining the problem considered by Morse and Kimball in 
Ref. 9* Their treatment is presented herein to facilitate the under- 
standing of Weiss’ development of the necessary and sufficient conditions 
contained in Appendix A. 

A. THE MORSE AED KIMBALL PROBLEM 

Morse and Kimball (Ref. 9 ) examined the problem of determining the 
firing pattern that maximizes the probability of at least one hit on a 
target when N rotmds are fired in a single salvo. An approximate 
solution is determined for large patterns. 

Morse and Kimball define the probability that a projectile will hit 
the small area element dxdy as f(x,y)dxdy. Then 

CO 00 

J J f(x,y)dxdy = N, (22) 

.00 .00 

where N is the number of rounds fired. The function f(x,y) is considered 
to be the pattern density fimction. Extensive changes in the hit 
probability can be obtained by changing the firing pattern. The expected 
number of lethal hits on a target located at coordinates x,y is L f(x,y), 
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where L is the lethal area of the target. Morse and Kimball approximate 
the expected nvimber of hits on the target using the Poisson probability- 
distribution. The probability of at least one hit is 1-exp-j^Lf (x,y )j-, (23) 

which is the probability of destroying the target. The total probability 
of destroying the target is then 






(24) 



where f^ is the probability density for aiming the pattern, usually the 
normal density. The problem then is to determine the fimction f(x,y) 
which maximizes P subject to equation (22). The problem is maximize 



^ = J J (^»y)<^<iy, 



Subject to J J f(x,y)dxdy = N, (25) 

^00 ^00 

f(x,y) ^ 0. 

It is the method of solution of this problem by Morse and Kimball that 
is interesting. They consider a pertubation of the ftinction f(x,y); that 
is, an increase of f(x,y) by a small amount 6 at the point x^,y^, and a 
decrease of f(x,y) by the same amount 6 at x^,y^. The constraint remains 
satisfied while the objective function changes by an amount 

’^1^fp(x^ ,y^) - e"^^^^2*^2^fp(x2,y2))s dxdy. 

Suppose f(x^,y^) and f(x2>y2) ^xe > 0; then if 

e"Lf(xi,yi)f^(x^,y^) > 

P can be increased by a 6 > 0. Conversely, if 



rLf(xi,yi)f (x^,y^) < e"^^^’^2^f (X2,y2), 
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then P can be increased by 6 < 0. Hence, for that function f(x,y) 
which maximizes P, it must be the case that 

for all points x,y where f(x,y) > 0, For all such points 



-Lf(x,y) 



f (x,y) = c > 0. 



P 



( 26 ) 



(27) 



Now, instead of x^,y 2 a point where f(x,y)=0 is considered. If 

f(x^»y^) is decreased by 6 (which now must be positive) and f(x^,y^) is 
increased by 6, then the increase in the objective function is 






(28) 



which equals (^f^(x^,y^) - cj 6 dxdy. This implies a positive increase in 
P is possible if f (x, ,y,) > c. Hence, f(x,y) cannot equal 0 unless 
f^(x,y) ^ c. It is seen that the solution is 



1 yf (x,y). 
i(x,j) . I Ln (-^-j ) 



= 0 



if fp(^»y) > c, 
if fp(^»y) ^ c. 



(29) 



The "uriknown constant c is determined from the constraint 

1 / 

" ^ T— ) ■ 

fpXx,y)>c 



(30) 



The set of conditions of equations (29) form a set of necessary and 
sufficient conditions for determining that pattern that maximizes the 
total probability of destroying the target, 

Weiss* analysis of the problem of minimizing the fraction of targets 
surviving a salvo of N rounds is presented in detail in Appendix A. The 
original paper (Ref, 7) is greatly condensed and the entire development is 
presented in this paper to facilitate the understanding of the techniques 
used in the solution. 
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rv. THE SOLUTIOII 



A. ANOTHER REPRESENTATION OP THE PROBLEI^ 



Weiss' coverage problem discussed in Appendix 1 can be written as 

00 CD 

Minimize $ = / J Pg,(x^»y^) e~^^^t ’^t^dx^ dy^; 

_00 .CD 



Subject to 



.00 .00 




dy^ = N 

7a) *0 



where 

00 00 
.00 .00 



( 31 ) 



P^(x, ,y.) and P.(x ,y ) are the target and aimpoint distributions 
JL u ”U A 3, 3* 

described earlier. It is noted that in addition to the constraints 
above, it is trivially true that 

S J* PTK’yt)=^‘ (52) 

.00 .00 

The solution to this problem is determined in Appendix 1 to be 
.00 .00 

> C if Pj4(^a-5'a)'° 

where q(x^-x^^ ,y^-y^^ ) relates to the lethality function. It is not 
possible to simply substitute other than trivial target distributions or 
lethality fimctions into the necessary and svifficient conditions of 
equations (33) aJid determine the distribution of aimpoint s that 
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minimizes the number of survivors. However, a solution to these 
equations may be obtained by using a finite difference approximation to 
the integrals. This paper assumes that the solution to the finite 
difference approximation converges tC the solution of the original 
problem. The proof of this convergence, however, is beyond the scope 
of this paper. It is realized that finite difference techniques may 
guarantee solution for specific values but usually preclude the seeing 
of relationships between equation parameters. 

Equation ( 51 ) can be approximated as 



k k _ r ') 

Minimize T. E P^(x^,y^)e”^''^t’^t^ , 



k k 



(5i) 



Subject to E E P.(x ,y ) = N, 

- ^ li. Qi Si 

k k 

where ii(x^,y^) = E^ P^(x^,yjLnq(x^-x^, y^-y^ ) 






and the upper limit of summation, k, is the number of increments in the 
approximation. Equations (53) similiarly be approximated by 



( 35 ) 

>C if 



The corresponding finite approximation to equation (32) is 
k k 

E E P^(x^,y^) = 1. (56) 

It is noted that equations (55) are the derivatives with respect to 
P.(x ,y ) of the objective function of equation (34) when the expression 

.n. 3, 3. 
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for u(x^,y^) is substituted into the objective function. Equations ( 55 ) 
can be vrritten as 









= c 
> c 






(37) 



It is noted that the objective function is convex in P.(x ,y ) and 

XX 3. 3/ 

that equations (35) (57) are simply the results of Gibbs Lemma 

as applied to equations (34). 

A solution to equations (35) can be obtained with a computer by 

dividing the target area into cells siofficiently small for approximation 

of the small areas dx.dy, and dx dy . To simplify the computer procedure, 

the following notations will be used; 

= The expected fraction of the target in cell i; it is 
analogous to the target distribution. 

y. = The proportion of the total rotinds aimed at cell j; it is 
analogous to the aimpoint distribution P.(x ,y ); 

XX 3 3 

b. . = The probability that a target in cell i is killed by a 
round aimed at cell j; it is analogous to the lethality 
function q(x^-x^,y^-y^) . 

The problem in equation (34) can then be rewritten as 

k k 

k -E E y b 

Minimize $ = Z P.e i=1 j=1 ‘J, 

i=1 ^ 



Subject to E P . = 1 , 
i=1 ^ 

k 

E y . = 1 
0=1 ^ 



(38) 



1 



P. , y . , b. . ^ 0 
1’ 0 10 



Gibbs Lemma is discussed in Ref. 10, 
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The solution is in terms of the y.’s and k is the total number of cells 

0 

Equations (35) and (3'7) can be rewritten as 




= C if y . > 0 
3 

> C if y. = 0 
0 



( 39 ) 



or, equivalently, 




Y 



if y > 0 
J 

if y = 0 
J 



(40) 



< Y 



where Y = -C. 

The problem thus presented is similar to one presented by Professor 
John M. Danskin to his Games of Strategy class at the U, S. Naval 
Postgraduate School in December, 1970. 

B. THE COMPUTER SOLUTION TECHNIQUE 

A solution technique is to divide an area of one Kilometer square 
into cells of twenty meters square. The P^*s are then determined for 
each of the 2,500 cells from the target distribution. To solve the 
problem of firing at targets obeying Koopmans* moving target distribution 
described earlier, values of time of impact and speed of the target are 
inputs and. the P^*s are determined from equation (2) for each cell of 
the simulated target area. The b. .’s are determined from the lethality 
fvinctions for the weapon system being investigated. For an artillery 
weapon it is the case that the ^ cells j more than, 

say, 60 meters from cell i. This means the number of machine calculations 
required for solving the problem in equations (40) would be reduced 
considerably as the sum over j would be limited to the 48 cells 
immediately surrounding cell i and, of covirse, cell i itself. 
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It is noted that if b. . =0 for all j ^ i; that is, the probability 
that targets outside the aimpoint cell are killed is zero, the solution 
is 



1 



P.b. . 



y Ln 

^ ij - Y 



if P.h. . > Y , 
1 10 ’ 



= 0 



if P.h. . Y . 
1 10 



Y is determined from the equation 



k . P.b.. 

E T-1- Ln — — ^ = 1 
i=1 ’’ij V 



which is analyzed in Ref, 10, 

After the values for the P^'s have been initialized and the lethality 

function formulas are determined for the cells surroTinding each cell i, 

a beginning set of values for each y. is initialized, A possible set of 

0 

j^'s is y.j=1 , y^=,,,= y2500~*^‘ then necessary to determine if that 

initialized set of y.'s give a minimum value of $ and, if not, in which 

0 

direction to move (i,e, what other vector of y.'s) to produce a lower §, 

0 

An algorithm has been developed by Professor Danskin that will cause 
convergence to a vector of y.'s that produce a minimum $, The computer 
program for this algorithm as v/ell as the application of the algorithm 
to the problem discussed in this paper are presented as a separate thesis 
for Professor Danskin by Major Paul T, Zmuida, USA, 
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APPENDIX A 



Weiss desires to minimize 
00 00 



« 0 O —00 2.=1 

Subject 'to ^ ^ 

J J P^(x^,y^)dx^dy^ = 1, 

—00 —00 

where $ is the expected fraction of the target sirrviving, P (x ,y.) is 

1 X 

the taxget distribution, and n q(x. -x . ,y . -y . ; is the probability a 

, . X ai u ai 
1=1 

target at coordinates survives a salvo of N rounds when the ith 

round is aimed at aimpoint x . ,y Weiss wishes to minimize $ by a 

ax ax 

proper choice of x . ,y . . 

ax ax 

A new function is defined such that 
N 

u(xt,yt) = -Ln _TT q(x^-x^., y^-y^p, 

X=| 

= -[in +...+ Ln 

= -.2. Ln q (^t-’^ai-^t-5'ai)’ 



i=1 



N 



or -u(r^,yp Ln q 

The objective fxinction can then be rewritten as 



(A-2) 



00 —CD 



(A-5) 
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A: 




Weiss then makes an approximation that a distribution function for the 



aimpoints can be substituted for the exact knowledge of x . ,y . so that 

3^3. 



the number of rounds aimed at area dx dy is P,(x ,y )dx dy , 

a a A'' a a"^ a a 

number of rounds is 
00 00 

J' / = ’*> 

mOO .00 

and since negative rounds cannot be fired, 

= 0 . 

Equation (A-2) can then be written as 



The total 



(A-4) 



(A-5) 






(A-6) 



The complete problem was minimize 

^ - S 5 e"'"K’^t)dx^,dy^, 

.00 .00 

subject to 
00 00 

X X = ■' 

.00 .00 

x"x 

.00 .00 

Pj,(xt,yt) s 0, ^ 0 

Weiss continues his solution using techniques similar to those of Morse 

and Kimball described earlier. For any arbitrary P.(x ,y ) add a small 

A a a 

increment 6 to P.(x ,y ) over the interval Ax at x and Ay at y , 

A a a a aL 

Now examine the change in $, The procedure can be illustrated as shown 
in Figure 5. 
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Figure 3* Illustration of small increment 6 

added to P.(x ,y ) over interval Ax at x ^ and Ay at y . 

ix 3/ 3/ 3/ QS. 3/ 3X 



The change in $ can be expressed as 

^ p /■ ^ = A § 

P.(x ,y ) +Ap (x ,y ) '^A^^a’^a'^ 

A' a’*' a'' A' a’’^a'^ 

From equation (A-5) , 



(A-7) 



J PiCxt-yp exp {- J J (PA(=^ 5 ,.yP+APA(=^ 5 ,.yp)Ln<l(x^-x y^,-y ) 

^00 ^00 _00 _00 



’dx 









A§ 



00 00 00 00 

^00 »00 ..00 .00 

00 00 , 00 00 

-J ^ {-I J" 



.00 .00 



= f f P„(=<*.yJ ^'A^'a-l'a'^ltVal'^t-l'aP'^a'^^a 

_ J. U "t _00 “ ® 

—00 _00 •" 

.[e- L L^A(=^a-^a)^ ^(=^t-al'^t-=^al)^aVl] ^^dy,, 
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1 

\ 






- J" J 



e 



T 



-' I '^t'^f 



(A-8) 



Since P^(x^,y^)=6 and Ln q (x-x^^,y-y ) is approximately constant over 



a very small area dx dy , 

a a 









-ijdx^dy^. (A-9) 



If the bracketed term is now expanded in a Taylor's expansion of e and 
a first order approximation is applied, then 

" “ J" L l(^t-^3i.y^-y3l)Jdx^dy^ 

»oo 

“ “ -u(x^,y^) 

_00 _CD 

Weiss defines 

® ® -u(x^,y^) 

-IS PT^^t’^t)® ^(^t-^al’^t-yal^'^t'^^t = P^^al’^al^’ 

..00 ..CO 

and then 



A$ =“ PC^al’^al^^^V^a* 



(A-10) 



Now either P.(x ,y ) >0 orP(x ,y ) =0. That is, either some roirnds 

±L Si Si A Si Si 

or no ro\mds are aimed at the immediate vicinity of x ,y . If 

Si Si 

P (x ,y ) >0, 6 can be added or subtracted from it; if P (x ,y ) =0, 

Si Si XX 3/ 



6 can only be added to it. 



Weiss considers two points ^hd ^g^2’^a2 ^A^^a’^a^ ^ 

6 is added to P,(x ,y ) at one point and subtracted from the other. The 

±i Si Si 

change in $ can be determined from equation (A-10}. 



24 



(A-1 1 ) 



A$ (x 



a1 






“ [ P (\l ’^'al)- f ^''^^a^a- 



If (A-11) is positive then $ can be reduced further by changing the sign 
of 6, which it is permissible to do since P.(x ,y ) > 0 at both points. 

.A. 3. 

It is clear that if 



P <==a1-ya1^ ^ 


P (^a2’^a2^ A$ < 0 


if 6 > 0, 




P ^^a2’^a2^ A$ < 0 


if 6 < 0. 



The conditions for an optimum exists v.'hen 

P ^^al’^al^ " P ^^a2’^a2^ " constant C, (A-12) 

for all points where P.(x ,y ) > 0, 

A 3. 3« 

Weiss then examines a point x ,,y ^ where P.(x ,y ) = 0, P.(x ,y ) 

3/^ Quj A 3/ 3. A 3. 3. 

can only be added to at this point. Comparing "to "the point 

where P^(x^»yJ > 0» 



A» = [p (\5.ya3) - P(^a1’^al0 **'a'^a' 

which is negative if and only if 

That is to say $ can be decreased only if p C. At the 

optimum solution § can no longer be decreased. 

The original problem is now reduced to solving the following conditions 



P Kl-^al) - 
< c 



“ ^tSV-rs} =■ ° 



(A-1 4) 
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I 



where »^a1 ^ " I S ^K"^a1 ’^t-^al 

»00 ^00 

and u(x^,y^) is as defined in equation (A-6). 

Weiss has determined necessary and sufficient conditions for the 
aimpoint distribution to be positive. The problem now has been reduced 
to solving equations (A-I4) which is a difficult problem and to which 
there exists no general solution in closed form. 
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